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Abstract 

The  analytical  sulutiuu  of  the  Mean  Spherical  Appruximatiuu  for  the  case  of  equal 
size  ious  aud  different  size  solvent  is  reexamined  using  only  two  parameters  :  a  polar¬ 
ization  parameter  A  and  a  screening  parameter  F.  We  show  that  the  ion  dipole  cross 
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energy  parameter,  which  iu  previous  work  was  obtained  solving  a  cubit  equation,  can 
be  obtained  from  a  linear  algebraic  equation.  Therefore,  the  inverse  problem  of  cal¬ 
culating  the  reduced  charge  parameter  do  und  the  reduced  dipole  parameter  d-i  from 
A  and  F  is  reduced  to  a  system  of  two  equations:  a  cubic  for  do,  and  a  linear  for  d^. 
Simpler  expressions  fur  the  thermodynamic  parameters  are  also  obtained. 
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1  Introduction 


Cousiderable  advances  have  been  made  over  tha  last  thirty  years  in  the  statistical  me¬ 
chanical  description  of  electrolyte  solutions  in  molecular  solvents.  This  work  has  been  carried 
out  using  different  approaches:  computer  simulations,^’^  accurate  integral  equations,  such  as 
the  hypernetted  chain  equation,^  which  were  applied  to  this  case  by  Patey  and  coworkers'*- ^ 
and  simple  analytical  theories,  based  on  the  mean  spherical  approximation  (MSA).*’  The  last 
approach  is  attractive  for  chemists  examining  the  thermod3Tiamic  properties  of  electrolyte 
solutions  not  only  because  the  model  gives  rather  simple  analytical  results  but  also  because 
it  shares  with  the  Debye  Huckel  theory  the  remarkable  simple  description  in  terms  of  a  single 
screening  parameter  for  any  arbitrary  mixture  of  electrolytes,  and  has  the  added  bonus  of 
satisfying  the  large  charge,  large  density  limits  of  Onsager.^’**  In  the  present  work  we  consider 
a  mixture  of  a  salt  of  equal  size  hard  sphere  ions,  and  a  solvent  which  is  represented  by  a 
hard  sphere  of  different  size  with  a  permanent  point  dipole.  Early  work  at  the  non-primitive 
level**’  ***  was  restricted  ,  for  technical  reasons,  to  the  case  of  ions  and  solvent  of  equal  diam¬ 
eter.  The  extension  of  the  analytical  result  to  mixtures  of  arbitrary  size  ions  and  solvents 
is  considerably  more  complex**"* '*  but  the  remarkable  fact  that  the  excess  ionic  properties 
depend  on  a  single  scaling,  Debye-like  parameter  is  still  retained  by  this  approximation.  The 
equations  for  the  most  general  case  appear  to  be  rather  complex,  but  the  semirestricted  case 
of  equal  size  ions  of  diameter  ai  and  a  different  size  dipole  of  diameter  a,  is  both  interesting 
and  tractable.  In  recent  years**  it  has  become  apparent  that  the  underlying  structure  of  the 
MSA  consists  of  a  number  of  scaling  parameters  which  is  equal  to  the  number  of  independent 
interaction  parameters  of  the  probleUi.  For  example,  for  the  ion-dipole  mixture  there  are 
only  two  interaction  parameters,  the  charge  and  the  dipole  moment,  and  correspondingly, 
there  will  be  oidy  two  scaling  parameters,  F,  for  the  charges,  and  A  for  the  '^ipoles.  For  the 
primitive  model  of  ionic  solutions  in  the  semirestricted  case**’the  parameter  F  is  determined 
from  the  equation 

twkuT 

where  tho  ionic  charge  is  z,t  and  number  density  p,  =  M,/V,  where  A/i  is  the  number  of  ions 
and  V  is  the  volume  of  the  system.  The  temperature  is  T  and  Boltzmann's  constant  is  ky. 
For  the  case  of  only  one 


?T 


-f-  F(7, 


=  4F' 


(1) 
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iou  diameter  cr,  the  screeuiug  parameter  F  is  related  to  the  Debye  screeuiiig  parameter 


K 


i 


iirt^ 


(2) 


by 


(1  +2ra)''  =  (1  +  2/ca) 


(3) 


For  a  system  of  spheres  of  diameter  a,  aud  a  permanent  dipole  moment  /i,  the  MSA  result 
tan  be  expressed  in  terms  of  a  single  parameter  A  .  Following  Wertheim,'^  we  have 


A'(A  +  2)'„ 

“i  —  n  (1 


(4) 


where 


4  = 


iTrp^n'i 


(5) 


3k  yT 

and  /),  is  the  solvent  number  density.  Furthermore,  the  MSA  dielectric  constant  e^v  of  the 
solvent  is  given  by 


ew  = 


A"(A  +  1)" 

16 


(6) 


As  has  been  often  done  in  the  literature,  the  parameter  A  can  be  computed  directly  from 
the  dielectric  constant  tw  using  the  above  cubic  equation.  This  parametrization  defines  an 
effective  polarization  parameter. 


This  paper  is  the  continuation  of  previous  work  of  Blum  aud  Wei^'^  aud  of  Blum 
aud  Fawcett.^^  The  notation  of  Blum  aud  Wei  will  be  used  throughout.  We  show  that  for 
the  semirestricted  case  of  the  ion-dipole  model  of  ionic  solutions  an  ion-dipole  interaction 
parameter  bi  is  obtained  from  a  linear  equation  whic  is  also  a  fuctiou  of  the  reduced  charge 
do-  In  our  previous  work^**’^^  this  parameter  was  obtained  as  a  solution  of  a  cubic  equation 
with  coefficients  that  depend  on  F,  the  dielectric  constant  tyy  aud  A.  When  cr,  0  we 
recover  the  primitive  model  MSA  equations,  aud  the  hydration  thermodynamics  is  that  of 
the  continuum  model  as  described  by  the  Born  equation.  In  the  next  section  we  discuss  the 
new  form  of  the  MSA  equations  for  the  semirestricted  case  of  equal  size  ions  aud  a  different 
size  hard  dipole.  The  thermodynamics  for  this  system  is  derived  in  the  final  section  using 
the  expressions  of  Vericat  et  al.'*  aud  of  Blum  et  al.*'* 
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2  Theory 


We  suiuiuatize  the  results  of  previous  work.'*’'*’''^  We  use  the  iuvariaut  expausiou 
formalism, in  which  the  total  pair  correlation  h{12)  is  expanded  in  terms  of  rotational 
invariants 

/i(12)  =  h'^{ry,)  +  (7) 

where  h”'’'*{ri2)  is  the  coefficient  of  the  invariant  expansion,  which  depend  only  on  the 
distance  ri2  between  spheres  1  and  2.  The  rotational  invariants  depend  oiJy  on  the 

mutual  orientations  of  the  molecules.  For  the  present  case  the  relevant  correlation  functions 
are 


•  ron-ion: 

/.„(.)  =  (1/2)  [A!i!?(r)  -  fc5^^(r)] 

(8) 

•  ion-dipole: 

h^r)  =  (1/2)  -  &’.'(>■)]  (f  •  M 

(9) 

•  dipole-dipole: 

hnn(r)  =  -v^ii‘„"(r))i,  -  it; 

+  yyAin  (>•)  [3(f  •  (ii)(f  •  it',)  -  ill  ■  ill] 

(10) 

where  /t  is  the  unit  vector 

of  the  ’energy’  parameters 

in  the  direction  of  fi.  The  solution  of  the  MSA  is 

given  in  terms 

•  ion-ion: 

too 

6o  =  27r/?i  /  drhii{r)TCTi 

Jo 

(11) 

•  ion-dipole: 

b\  =  drhi„{r)(T,a, 

(12) 

•  dipole-dipole: 

=  dr  y„i 

(13) 
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which,  as  will  be  shown  below  are  proportional  to  the  ion-ion,  ion-dipole  and  dipole-dipole 
excess  internal  energy. **  In  the  MSA  they  are  functions  of  the  ion  charge  and  the  solvent 
dipole  moment,  through  the  parameters 


“o  "  koT  *  ^ 

(14) 

and 

J2  47r/9,/ii 

"  3kiiT 

(15) 

These  parameters  are 

required  to  satisfy  the  following  equations^^ 

ai  +  al  =  do 

(16) 

1 

a 

1 

3 

a 

II 

"a 

(17) 

+  [1  ~  —  yi  +  dit 

(18) 

Since 

I<m  =  {<rJa,)^[l  +  aiA] 

(19) 

and 

1  -  Kii  =  -  a2biA{at/2cri)] 

(20) 

with 

A  =  (1/2)(1  +  bv)  +  0ii{cr,IZa,) 

(21) 

can  be  rearranged  to 

«1  + 

(22) 

-  a-i^i  =  doA4 

2cTj 

(23) 

^  +0i  =  -H  A"^" 

(24) 

where  we  have  used  the  definitions 


«.  =  1 

(25) 

(26) 

Dy  -  2  [^,(1  +  4»)  ] 

(27) 
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(28) 


a.ud 


which  cau  also  be  writteu^^ 


A=d2 


hid{}G  f 
4A(7i 


1  +  60  + 


jdacr. 


ScTj 


A  —  ^2  —  doA; 


h\<T, 


w 


here 


with 


20iVa, 

D  =  1  +  Si 


„  _  4;  4;{A  +  2)' 

°'-wr 


Also 


=  1  -  6.2/6 
/?3  =  1  +  6.2/3 

=  1  +  6.2/12 

It  will  turn  out  couveuieut  to  use  the  parameter 

/33 


Theu 


A  = 

6.2  =  6 
/?3  =  3 

y9«  = 


/?« 

A  -  1 
A  +  2 
A 

A  +  2 
3 


A  +  2 


_  _  A(A  +  2) 

3 

t4(A  +  2) 
3{A  + 1)’^ 

Futhermore  we  have  the  dielectric  coustaut*’^ 

yi  + 1)^ 


€W  =  —  = 

vt 
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(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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As  has  beeu  previously  showu*^  *^  Eq.(23)  cau  be  replaced  by 


UibiCTt 


+  =  doAyi 


which  together  with  Eqs.(22)  aud  (24)  form  a  system  of  three  equations  for  the  three  un¬ 
knowns  6i),  b\,b-i.  However  these  equations  are  still  complex.  A  much  simpler  set  of  equations 
is  obtained  when  we  use  the  proper  scaling  lengths.'*’’'^  We  define  F  through  the  relation 

t  ,  h{  1  ,  /3bO-,1  _  -r<7.  ^  2BiK 

^  l-|-r<Ti  4^^  l  +  T(Ti  ZcT,  1  +  Tai  V 


ra.  =  - 


(1  +  b,-?!^) 


and  A  from 


0ti 


After  some  long  but  straightforward  algebra  we  get 

^  n  /t  I  ■p  \  /^sFcTi 

=  pr<7,(l  +  r<rO  =  — 


Ui  =  -ai 


6i  _  biTcr, 
2^  ~  ~2Dj.  O 


where  now 


2(1  -I-  Fa.) 

where  we  have  defined  the  convenient  parameter 


where 


1  +  Fa,  -H  ^  1  +  Fa,  -H  f 

Cl 


f  =  — 

Aa, 


A  =  (1/2)[i  +  ^  +  ^]=(P/2)[j^+^ 


36 


Furthermore  we  get  for  Eq.(29) 


4  J  J  /^*\/  / ^*  \ 
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Also 


1  -  Kn  =  (yo/P) 


<^0-77  [c/  + 

Si 


1  +  5i(6/0)r(7,  i  +  (i  +  rcr,) 


.3iiCr, 

3cr, 


With  this  the  boundary  condition 

25i  r  ,, 


the  MSA  equation  (43  )reduces  to 


bidodi  —  d^{a,lcr,)h.  +  (A/^,)(l  +  Ci6)]  • 


2r<T.(i  +  r»,)(i-(i>-i)({,/C.)l  = 


y'-W 


together  with  the  cubic  equation  for  the  ion-dipole  parameter  i>i 

These  can  be  written  in  a  more  convenient  form  using  Eq.(59) 

do(l  +  Si)^ 


2r<T.(i-hrai)(i~Si{6/c.)]  = 


Eq.(60)  is  a  cubic  equation  in  Bi 

(1 +B, )'[('  + 0,]  =t,»  [C.-B,C1' 


(55) 

(56) 

(57) 

(58) 


(59) 


(60) 


(61) 


(62) 


Eliminating  (1  between  these  two  equations,  we  get  a  linear  expression  in  Bi  (  in  place 

of  the  cubic'^),  which  is  easily  solved  to  yield 


01  =  c 


[1  - 


(63) 


I'i  + 

This  is  one  of  our  main  results;  The  parameter  61  is  now  explicitly  obtained  as  a  function  of 
the  other  parameters  of  the  problem,  where 

2r<7,(l  -t-  Fa,) 


I/,  = 


Furthermore 


1  +  Si  = 


duy/fW 

I'.  +  6C. 


(64) 


(65) 
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1  -  (6/C.)Si  = 


cm+g.) 

Vi  +  ^iCi 

Vi{l  +  ^.0) 


Vi  +  (,Q 


=  [1  +  B^Y  =  twv] 


,(i  +  6C.) 


V,  +  (f.Ci 


so  that  from  eq.(61)  we  get 


2r(7.(i  +  =  do!^ 


■V,  -  ^iC,  -  Ct^  +  viCY 
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Vi  +  (,C, 


Substituting 


We  have  that 


Vi  +  6C. 

vtev^il  +  ^lOYvi  +  6C.)  =  {-V.  -  6C.  -  C"  +  v,(fY 

(1  +  6)(i  +  rai)rcr, 


Vi  +  60  = 


1  +  00  = 


T(Ti 

1  +  r<Ti  + 1/0 


1  +  Ta,  + 1/ 0 


0  +  2 


(1  +  r<Ti)(l  +  Fcr,  +  1/0) 


do\/tw 


+  00  +  0  ''.o  (i  +  r<T.  +  i/o)" 

Finally,  after  some  leng^j^  algebra 


1  + 


2(l  +  2ra. +  1/0) 


<fov/«wO 


do(l  +  0)(2  +  20  +  4rcr,0  +  doy/tw^i ) 

=  4r<T,Vcw'(i  +  rcr,o)(i  +  0  +  r^^iO)^ 

(2  +  2r<7i  +  20  +  4r<r,0  +  2rV/0  +  di,yjtwii) 

which  can  be  written  in  the  form 

[do  +  2(1  +  1/0  +  2r(7,)/(ye^v0)'" 


do 


(1  + 1/0  +  r<^») 


=  +  2(1  +  r<7.)(i  + 1/0  +  r<7.)/(0ve^v)] 


This  is  a  cubic  equation  for  the  parameter  do- 


^66) 

(67) 

(68) 

(69) 


(70) 

(71) 

(72) 


(73) 


(74) 


do[do  +  Cj]^  —  c-i{do  +  ca) 


(75) 


with 
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_  or<7.(i  +  rcr.) 
C«6\/^i^ 

A  second  equation  is  laboriously  obtained  tom  Eq.(18) 


(76^ 

(”) 

(78) 


1  +  (P/2)  Ad. 


1  0tcr» 


1  +  Fcr,  3ai 


+  :r7: — Ad, - ,  .  ,  .  t, - =  Vi  +  dj 


where 


Ad.  = 


1  +  1/^i  +  FcTi 
2r<7,(l  +  Fcr.)  do 


_ 


v^ew 


(79) 


(80) 


(81) 


From  these  two  equations  the  inverse  problem  can  be  solved  explicitly,  that  is,  to  compute 
the  parameters  do  and  d-i  from  F  and  A.  When  the  solvent  diameter  cr,  goes  to  zero  we 
recover  the  correct  primitive  model  result 

do 


2Fa.(l  +  Fa.)  = 


\/Fw 


(82) 


which  is  equal  to  the  primitive  model  MSA  Eq.(3)  in  the  form 

2Fa(l  +  Fa)  =  KO 


(83) 


with  the  MSA  Wertheim  dielectric  constant  tw 


3  Thermodynamics 

The  excess  internal  energy  for  the  semirestricted  case  of  equal  size  ions  and  a  different 
size  solvent  is  similar  to  the  completely  resticted  case  of  all  equal  size  spheres.  The  thermo¬ 
dynamic  relations  for  the  completely  restricted  case  were  first  derived  by  Vericat  and  Blum.** 
The  generalization  of  the  thermod3atamic  expressions  to  the  completely  general  case  and  the 
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semiresticted  case  were  obtained  later  by  Wei  and  We  use  the  reduced  quantities  of 

the  second  work.  The  internal  energy  is 

afEliVkuT)  =  ^  d^6o  -  2  dod,6i  -  2  djh-.  (84) 

the  Helmholtz  free  energy  is  obtained  using  the  method  of  Hoye  and  Stell'^^ 

oiAHykaT)  =  (2dX  -  2dudA  -  J')  (85) 

d'  =  m‘  +  [l  +  (^)]  (~)  IQUf  +  (^y  [[Q'J‘  +  2(9')']  (86) 

It  will  be  convenient  to  write  the  Helmholtz  free  energy  as  the  sum  of  a  internal  energy 
and  an  entropy  term 

AjiVktiT)  =  E/{VksT)  -  S/{Vkti)  (87) 

where 


-  afS/iVk^)  =  -dS6o  +  4(^)dod,6,+6(^) 

Using  the  tesults  of  Wei  nud  Blum/^  Eqs.(16-20)  nud 

,  _  1  -  82/24  _  „(A-l)(A+3) 

’  '(1+V12)»  (A  +  1)' 


Q',  =  -  2  +  A/flf  =  -|irv,'  +  e.i 

QU  =  ^(1  +  r»,)(i  +  r9,«,) 

=  2  [(l/A)  (/Jj  -  (£-)  6.<i.(3A  -  2I>, ))  -  l] 
After  some  leughty  calculations  we  get 

-  S/(Vk^)  =  ^  [r"  +  ^  ([1  +  -  3)]  +  w 

OTT  [  (T;^  \  Cw  / 

where  the  coefficients  T'  is  a  second  order  polynomial  in  the  variable 

_ M _ 

V‘‘  2(i  +  r<7.)(i  +  (n-rtT.)^i) 
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where  M,  has  been  defined  in  eq.(80).  When  this  variable  is  zero  then  the  system  behaves 
like  a  sum  of  a  primitive  model  electrolyte  and  a  dipolar  fluid.  Our  expressions  reproduce 
the  known  results^^^®  for  the  ionic  part.  For  the  dipolar  part  we  get  a  new  and  simple 
expression,  which,  however,  is  in  full  agreement  with  the  result  of  Rushbrooke  et  al.^*’ 


We  have 

r  =  {\l\2Tr)[Tt,  +  Tj\  (95) 

with 

Tb  =  Tfco  -h  Tbi  +  Tb2  (96) 

Tm  =  -  [^]  4(r<T,)»(i  +  r<T4  2  - 1/('  +  B.  +  (i/r».)(i  +  b./c?)  (i  +  (i  + 


Til  =  16  [^]  r<i,(i  +  r<7i)  1  +  ^  +  r<’i(i  +  s./c?)  (i  +  (i  + 

(97) 

(98) 

Tb2 

=  (^)  66-4(i^  +  y't  -  y^] 

-“■©■(SI 

(99) 

where 

n  =  -1  +  «./c.)r.7i  (i  +  (1  + 

Furthermore 

Tj  =  Tjo  +  Tj\  +  Tj2 

(100) 

T  l[®'' 

Tj,  =  4 

(i-rV)[2rv,^  +  Si(i  +  rv.‘^)] 

(101) 

(102) 

Tn  =  (^y  [iQiJ'  -  4(A'  -  !)■'] 

=  -4  [^]  [a'  -  (r.T,/C,)'')  [(A'  -  1)(1  +  P)  -  B,  (l  -  (r<r,/C,)')]  (103) 

The  excess  pressure  can  also  be  computed.^*  The  expression  is*'* 


P/kyT  =  S/VkaT 


(104) 


G  =  E 


(105) 


Then, 

still  holds. 

The  thermodynamic  expressions  for  the  infinite  dilution  limit”  was  discussed  by  Chau 
and  colaborators,^^  and  in  a  more  comprehensive  way  by  Garisto  et  al.^**  Our  expressions 
agree  with  these  works  in  that  limit.  A  more  detailed  discussion  of  the  applications  of  these 
results  will  be  published  in  future  work. 
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